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1 Introduction 



The continuum theories of the flow of vehicular traffic are analogous to the 

hydrodynamic description of fluid flow. The kinetic theories of vehicular traffic [Q- 
pl are extensions of the kinetic theory of gases whereas the " car- following models" 
|j [To] and the recent particle-hopping models [pd|-p9[ are analogues of the models 
of driven system of interacting particles p0[ |. Therefore, the theoretical techniques 
of fluid dynamics and non-equilibrium statistical mechanics have turned out to 
be powerful tools in the study of a wide variety of problems in the science and 
engineering of vehicular traffic |u]-|33| . One of the most popular " particle- hopping 
models" of vehicular traffic is the Nagel-Schreckenberg (NaSch) model pH| , which 
may be regarded as a stochastic cellular automaton (CA) p4| , where the time- 
evolution of the vehicles is formulated in terms of a set of "rules" . Several extensions 
of these rules have been suggested in recent years (see, e.g., ||l]-|2g]) to make the 
model more realistic; in this paper we study some of the statistical properties of 
one such extension [25] . 

The flux versus density of the vehicles is known as the fundamental diagram |3l| . 
The distance-headway (DH) is defined as the distance from a selected point on the 
lead vehicle (LV) to the same point on the following vehicle (FV) ; usually the front 
edges or bumpers are selected ]|f| . The time-headway (TH) is defined as the time 
interval between the departures (or arrivals) of two successive vehicles recorded by a 
detector placed at a fixed position on the highway . The fundamental diagram, 
the DH and TH distributions, the distributions of the sizes of the traffic jams as 
well as the distribution of the gaps between successive jams in the NaSch model 
on idealized single-lane highways have been calculated Jll - 20 1 . Very recently, the 
update "rules" of the NaSch model have been extended by adding a step where 
the so-called "slow-to-start" rules (which will be explained in the next section) [ p5| - 
p8| are implemented. The fundamental diagram of some of these CA models with 
" slow-to-start" rules have been investigated (2^] . The main aim of this paper is to 
calculate, by a combination of analytical and numerical methods, the distributions 
of DH, TH, jam sizes and distances between jams in the CA model with the "slow- 
to-start" rule suggested by Benjamin, Johnson and Hui (BJH) |^5| to study the 
effects of the slow-to-start rule on the spatio-temporal organization of traffic. We 
interpret the novel features of these distributions in terms of a phase separation of 
the traffic in the BJH model. 

The NaSch model and the BJH model |25|], which was formulated by BJH 
by incorporating one specific type of "slow-to-start" rule in the updating scheme 
of the NaSch model, are stated in section 2 for the sake of completeness. The 
theoretical techniques, which will be used in our analytical calculations, and the 
numerical methods, which will be followed in our computer simulations, are also 
explained briefly in the same section. The probability distributions, namely the 
distributions of DH, TH, jam sizes and gaps between jams, with which we are 
mainly concerned in this paper, are defined in section 2. Our results for the DH 
distributions are presented in section 3. Our analytical as well as numerical results 
for the distributions of jam sizes and jam gaps are presented in sections 4 and 5, 
respectively, while those for the TH distributions are given in section 6. We make 
a quantitative estimation of the regime of validity of our analytical expressions by 
computing an appropriate correlation function in section 7. A summary of our 
results and conclusions drawn from these are given in section 8. Three appendices 
contain details of our analytical calculations. 
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2 Models and theoretical techniques 



In the particle-hopping models a lane is represented by a one-dimensional lattice of 
L sites. Each of the lattice sites can be either empty or occupied by at most one 
"vehicle". If periodic boundary condition is imposed, the density c of the vehicles 
is N/L where N(< L) is the total number of vehicles. Throughout this paper we 
shall follow the convention that the vehicles move from left to right (i.e., along the 
positive X-axis), so that in a given configuration (a, b), of a pair of sites, b refers to 
the LV and a refers to the FV. 

2.1 The models 

2.1.1 The NaSch model 

For the sake of completeness we briefly recall the definition of the NaSch model 
In this model, the speed V of each vehicle can take one of the V max + 1 
allowed integer values V = 0, 1, ■■.,V max . Suppose, X n and V n denote the position 
and speed, respectively, of the n-th vehicle. Then, d n = X n+ i — X n , is the gap in 
between the n-th vehicle and the vehicle in front of it at time t. At each discrete 
time step t —> t + 1 , the arrangement of N vehicles is updated in parallel according 
to the following " rules" : 

Step 1: Acceleration. V n — > min(V n + 1, V max ). 
Step 2: Deceleration (due to other vehicles). V n — > min(V n , d n — 1). 
Step 3: Randomization. V n max(V n — 1,0) with probability p. 
Step 4 : Vehicle movement. X n — » X n + V n . 

Non- vanishing braking probability p is essential for a realistic modeling of traffic 
flow Jl3| and, therefore, the NaSch model may be regarded as stochastic cellular 
automata 

2.1.2 The BJH model 

To our knowledge, so far three different versions of the " slow-to-start" rule have 
been formulated p5|-p8t ; in this paper we shall consider only the slow-to-start rule 
introduced by Benjamin et al. (BJH) p5| ). BJH modified the updating rules of the 
NaSch model by introducing an extra step where their " slow-to-start" rule is imple- 
mented. According this "slow-to-start" rule, the vehicles which had to brake due 
to the next vehicle ahead will move on the next opportunity only with probability 
1 — p s . The steps of the update rules can be stated as follows: 

Step 1: Acceleration. V n — * min(V n + l,V max ). 

Step 2: Slow-to-start rule: If flag = 1, then V n — > with probability p s . 
Step 3: Blockage (due to other vehicles). V n — * min{V n ,d n — 1) and, then, 
flag = 1 if V n = 0, else flag = 0. 

Step 4 : Randomization. V n max(V n — 1,0) with probability p. 
Step 5: Vehicle movement. X n — > X n + V n . 

Here flag is a label distinguishing vehicles which have to obey the slow-to-start rule 
(flag — 1) from those which do not have to (flag = 0). 

For p s = the above rules reduce to those of the NaSch model. Note that, 
in both these models, if V max — 1, a vehicle can come to an instantaneous halt, 
even at vanishingly low densities, because of random braking. In contrast, for all 
V max > 1, vehicles do not stop spontaneously at sufficiently low densities for which 
interactions among the vehicles is negligibly small. Therefore, effectively free flow 
of traffic takes place when the density of vehicles is sufficiently low whereas high 
density leads to congestion and traffic jams. The density, c , corresponding to the 
maximum flux is usually called the optimum density. At sufficiently high densities, 



3 



the B JH model with V max > 1 is known to exhibit phase separation where one of the 
two phases consists of a macroscopic jam while the other corresponds to free flow 
of traffic. Moreover, the particle-hole symmetry of the NaSch model for V max = 1 
is lost in the corresponding B JH model because of the slow-to-start rule |^8j ■ 

2.2 Theoretical techniques 

We compute the DH distribution in the BJH model with V max = 1 using some 
analytical results derived earlier psfl within the framework of a car-oriented mean- 
field (COMF) theory ]l4|] . Our analytical calculations of the distributions of jam 
sizes in the BJH model with V ma x = 1 are also carried out within the framework of 
the same COMF theory. We derive analytical expressions for the distributions of 
the gaps between successive jams in the BJH model for V max = 1 using a 2-cluster 
site-oriented mean-field (SOMF) theory which leads to the exact results for the 
NaSch model with V max = 1 in the limit p s — > 0, although the results for the BJH 
model (p s ^ 0) are not exact. We also present an alternative derivation of the 
distributions of the gaps between successive jams in terms of a hybrid approach 
which will be explained in section 5. Finally, we derive the analytical expression 
for the TH distributions in the BJH model with V max = 1 following an extension of 
the approach used earlier jl9) for the corresponding calculation for NaSch model. 

In our computer simulations we begin with a random initial configuration and 
let a system of L — 1000 sites evolve for 10 4 steps which is long enough to ensure 
that the system reaches the steady state; then we let the system evolve for 10 5 
further steps during which we compute the properties of our interest. Then the 
data are averaged over 100 samples (i.e., 100 initial configurations) for each of the 
sets of values of the parameters. We have analyzed data for a wide range of values 
of the parameters p and p s . However, we present here data mainly for the regime 
p <C p s where effects of the slow-to-start rule are the strongest. 

2.3 The distributions of our interest 

In this section we define the distributions of our interest. The number of empty 
lattice sites in front of a vehicle is taken to be a measure of the corresponding DH. 
Suppose Vdh(k) denotes the probability of a DH k; more precisely, Vdh(k) is the 
conditional probability of finding a string of k empty sites in front of a site which 
is given to be occupied by a vehicle. Note that Vdh(0) is also the total fraction of 
vehicles which are simultaneously in the jammed state. 

A jam of length k is defined as a string of k successive stopped vehicles, i.e. we 
are considering only compact jams. Similarly, when there arc k lattice sites between 
two successive jams, each occupied by a moving vehicle or is vacant then we say 
that there is a gap of length k between the two successive jams; we are interested 
in the distribution Vj g {k). 

We define the TH as follows fis|] . A detector is placed at one of the lattice 
sites. In order to measure the TH between a given pair of vehicles the detector 
is set to t = instantaneously when the LV of the pair leaves the detector site. 
From then onwards the detector counts the number of discrete time steps up to 
t = t when the FV leaves the same detector site; r is the TH between the pair 
of vehicles under consideration. As soon as the FV leaves the detector site the 
detector instantaneously resets again to t = in order to measure the TH between 
this vehicle and the vehicle following it. We denote the TH distribution by the 
symbol P t h(r). 
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3 The DH distribution 



Suppose, V(k) denotes the probability to find exactly k empty cells in front of a 
vehicle. Only vehicles which have exactly one empty cell in front might be affected 
by the slow-to-start rule p8| . We denote the density of these vehicles by 'P(l). The 
normalization requires 

X>(*)+£(i) = i. (i) 

fe>0 

Note that V dh {\) = V(l)+P(l), and V dh {k) = T(k) for all k + 1. The symbol g(t) 
{g(t)) is the probability that the vehicle moves (does not move) in the next time 
step; therefore, 

g = q^2V(k)+qq s V(l) (2) 

fc>l 

where q = 1 — p and q s = 1 — p s . Making the ansatz 

V(k) = Mz k (Jfc > 2) (3) 
one gets the following results ps|| : 

m= ^m_ qs P (1) (4) 



- _ qn o)(l-V(Q)) 

V{1) - l+ Ps qV(0) (5) 



where 



Af=i-^.[1-V(0)-V(l)-V(l)}, (6) 



Qi 1 - 9) 
q(l-V(0)) 

9 l+p s qV(0) (8) 
and V(0) is given as the root in the interval (0, 1) of the cubic equation 

cp 2 s q 2 V(0) 3 + q[qp 2 s (l-2c)+p s (l + c) + c}V(0) 2 

+ [qp s (l-3c)-2qc+l]V(0)-pc = 0. (9) 

In the special limit p a — > 0, 



c \ c 

qy 

g ->■ — ; 

c 

and 

z _^ (py/c) = 1 _ y 
1 - (qy/c) i - c 
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where 



V = 



1 



1 - y/l - Aqc(l - c) 



is the solution of the quadratic equation 

qy 2 - y + c(i - c) = o. 



(10) 



(ii) 



Moreover, in the special case p s = , the cubic equation (j|) for Vdh(0) reduces to 
the quadratic equation 

qcP dh {0 f + (1 - 2qc)T dh {0) - pc = 0, 

whose physically relevant solution 



„ / n x (2 g c-l) + v /(2 g c-l) 2 +4pgc T y 

= = 1 

2qc c 

is identical with the corresponding exact result [jl4|,[l7| for the NaSch model. Thus, 
in the limit p s — > 0, 



V dh (l)=V(l) + V(l) 



c(l - c) 



and 



c(l - c) 



1 - c 



fc-i 



for fc > 1, 



which are the corresponding known results for the NaSch model 0,117]] 
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Figure 1: The DH distributions in the BJH model with V m ax = 1 are shown on a 
semi-log plot for three different values of the density of vehicles c. The lines are 
obtained from the approximate analytical expressions for Vdh{k) while the discrete 
data points have been obtained from computer simulation. The common parameters 
are p = 0.05, p s = 0.50. 

For V m ax = 1 j and for given c, we have solved the cubic equation (^|) numerically 
for V(Q) and, hence, computed g, z,V(l),V(l),N, V{k) for all k > 2. Then, we 
calculated the distance-headway Vdh(k) for all k in the BJH model corresponding 
to the density c. A comparison of these results with those obtained from computer 
simulation for the same set of parameters (fig. establishes that for low densities, 
e.g. c = 0.10, and c = 0.25, the approximate analytical expressions for Vdh(k) 
derived in the COMF approximation, arc in good agreement with the simulation 
data. Larger deviations of the analytical results are observable for c = 0.5. 
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Figure 2: The DH distributions in the BJH model and the NaSch model, both with 
V max = 1 and p = 0.05, are shown on a semi-log plot for three different values of the 
density of vehicles c. The lines are obtained from the exact analytical expressions 
for Vdh(k) in the NaSch model while the discrete data points have been obtained 
from computer simulation of the BJH model with p s = 0.50. 




We have made a comparison of the DH distributions in the NaSch model and 
the BJH model, both with V max = 1 and with p — 0.05, in fig. ||. At the low density 
of c = 0.1 there is practically no difference in the DH distributions in these two 
models when p — 0.05. However, at the moderately high density of c = 0.25, large 
DH are slightly less probable in the BJH than in the NaSch model. Surprisingly, the 
difference in the DH distributions in these two models is less significant at c = 0.5 
than at c = 0.25. This, we believe, is a consequence of the fact that at a high density 
of c = 0.5 the DH distribution is decided mainly by the blockage of one vehicle by 
another, rather than the slow-to-start rule. Some typical DH distributions in the 
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Figure 3: Same as fig.|], except that V max = 5 and the lines are merely guides to 
the eye. 

BJH model with V m ax — 5 are shown in fig. ||. The common feature of the DH 
distributions in figs. and || is that for c — 0.1, which corresponds to free-flowing 
traffic, DH smaller than V max are strongly suppressed, while beyond V max the DH 
distribution decays exponentially. The density c = 0.5 is sufficiently high so that 
there is a very high probability for a vehicle to be stuck in a jam; this is consistent 
with the maximum of Pdh(k) located at k = 0. Similar peaks at k = have also 
been observed earlier [|l7| in the DH distribution in the NaSch model. What is more 
interesting is that the DH distribution in the BJH model with V max — 5 exhibits 
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peaks also at k = 5 and k = 10, i.e., at V max , and 2V max (albeit of gradually smaller 
heights) . This phenomenon can be understood by considering p = and recalling 
that, for the values of the parameters used in our simulation, phase separation of 
traffic into a " free- flowing" region and a region of jammed vehicles occurs at c = 0.5. 
If p = 0, then, the distances between the vehicles in the free flowing region of traffic 
are determined by the waiting times of the leading vehicle in a jam. If the leading 
vehicles leaves the jam in the very first opportunity then its DH is likely to be V max 
whereas its likely DH is 2V max if it starts at the next time step because of the 
" slow-to-start" rule. This perfect arrangement of the vehicles would correspond to 
delta function- like peaks at V max and 2V max , etc. but gets smeared out to sharp 
peaks, due to the fluctuations, if p is nonzero but not too large. 



4 The distribution of jam-sizes 



Within the framework of the COMF theory, the distribution of the sizes of the jams 
in the BJH model with V™, 



1 can be written as fl2C 



v js (k) = a/;- s [i - vmivm^ii - no)} 

where Mj S is a normalization coefficient. Since, the condition ^2u'Pj s (k) = 1 of 
normalization leads to Nj S = [1 — ^(O)] -1 , the normalized distribution of the jam 
sizes is given by 



p JS (k) = [i-v(0W(0W 



fe-i 



(12) 



In the limit p s 



0. 



P JS (k) 



y -(l- y - 



k-1 



c \ c, 

which is identical to the jam-size distribution in the NaSch model with V ma x 
where y is given by the equation (|To|). 
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Figure 4: The distribution of jam-sizes in the BJH model with V max = 1 are shown 
on a semi-log plot for the same values of c, p and p s as in fig. The lines 

are obtained from the approximate analytical expression ( |l2| ) for Vj S (k) while the 
discrete data points have been obtained from computer simulation. 



The jam-size distribution for the BJH model is well approximated by the ex- 
pression ( |l2|) at low densities of the vehicles as is evident from the fig. [| where 
comparison has been made with the corresponding numerical data obtained from 
computer simulation of the BJH model with the same set of parameters. How- 
ever, the approximate expression ( [l2] ) deviates more and more from the computer 
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Figure 5: The distributions of jam-sizes in the BJH model and the NaSch model, 
both with V m ax = 1 and p = 0.05, are shown on a semi-log plot for three different 
values of the density of vehicles c. The lines are obtained from the exact analytical 
expressions for Vj S (k) in the NaSch model while the discrete data points have been 
obtained from computer simulation of the BJH model with p s = 0.50. 



simulation data with increasing density for the densities we took into account. A 
comparison between the jam-size distributions in the NaSch and BJH models, both 
for V m ax = 1 and p = 0.05, is made in fig. ||. This comparison establishes that the 
difference in the jam-size distributions in these two models becomes wider with the 
increase of vehicle density c. At a high density of c = 0.5, for example, long jams 
are much more probable in the BJH model than in the NaSch model. A few typical 
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Figure 6: Same as fig. |J, except that V max = 5 and the lines are merely guides to 
the eye. 

distributions Vj S (k) for V max = 5 are shown in fig. ^. From figs. |, | and | we 
conclude that, for both V max = 1 and V ma x — 5, the jam size distributions Vj S {k) 
decay exponentially with the jam-size k. 
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5 The distribution of gaps between successive jams 
5.1 A hybrid approach 

Let us use a three-state variable s to describe the state of a site; s = — 1 if the site 
is empty, s = if the site is occupied by a vehicle with instantaneous speed V = 
and s = 1 if the site is occupied by a vehicle with instantaneous speed V = 1 |17j . 
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Note that this variable s does not distinguish between the slow-to-start vehicles and 
other vehicles which have instantaneous speed V = 0. We denote the probability of 
a configuration immediately after the slow-to-start stage by the symbol P(a, b ■ ■ ■ ) 
while that immediately after the randomization stage by p(c, d ■ • ■ ) and use the 
sequence 5 — 1 — 2 — 3 — 4 of the steps of the update rules, instead of 1 — 2 — 3 — 4 — 5. 

In terms of the state variable s, a 2 — cluster — like approximate decomposition 
of the probability distribution of the gaps between successive jams leads to the 
expression 



V jg {k)=N ji 



P 2 (Q | Sl)p 2 (si | S 2 ) " ■ -p 2 ( sfc-l | Sk)p2(sk_ | 0) 

{Si=±l} 



for k > 1 



(13) 



where Afj g is a normalization coefficient and the conditional probabilities p 2 (a|6) 
are defined through the equation 

i in P2(a,b) 

p 2 (a,0) +P2(a, -1) +P2(a, 1) 

Let us now define the 2x2 transfer matrix T with the following elements: 

T[l, 1] = p 2 (l | 1), T[l,2] =P2{1 | -1), 

r[2,l]=p a ( = l|l) ) ? 1 [2,2]=p 2 (^l|-l). 



In terms of this transfer matrix, the expression (|13|) for Pj g (k) can be recast as 

V jg (k) = AT is [p 2 (0 | l)T k -% 1] P2 (1 | 0) +p 2 (0 | -ljT*- 1 ^ L> 2 (1 | 0) 

+ P 2 (0 | l)!*"^, 2]p 2 ( = l | 0) + P2 (0 | -ljT*- 1 ^^^ | 0)] (16) 

where T k ~ 1 [i, j] refers to the [i,j] element of the matrix T k ~ 1 . It is straightforward 
to show that 

T*- l [l,l] = (A 1 A^- 1 -A 2 A^ 1 )/(A 1 -A 2 ), (17) 
T h ~ l [l,2] = (Xt 1 ~ A^/CAr - A 2 ), (18) 
T fe_1 [2, 1] = (Ax A| - X 2 X k )/(Ai - A 2 ), (19) 
T fe - 1 [2,2] = (Aj - A§)/(A! - A 2 ), (20) 

where Ai and A 2 are the two eigenvalues of the transfer matrix T. 

As shown in the appendix A, the 2-cluster-like conditional probabilities p 2 (a|6) 
can be related to the quantities like g, z,P(0),P(l),P(l), etc. which occur in the 
COMF theory. Using these relations we get 

T[l,l] = p 2 (l|l)=Q ! (21) 
T[l,2] = p2(l|-l) = l, (22) 

Troil , , m cg(l -P(0)) 

T[2,l] = P2( = l\l) = c[1 _ m]+z{1 _ cy (23) 

T[2,2] = P2 (-l|-l)= ^fi~ C li V ( 24 ) 

c[l - 7 , (0)J + z(l - c) 

and, hence @, with D = c[p + gP(0) + gp s P(l)], 

eP(0)g X^-Xt 1 c(l - g)[l -V(0)\ 



V jg (k) = Afjg 



D Ai-A 2 c[l-P(0)]+2(l-c) 

V k \k 



( cP(0)\ X k - X k c(l-, 9 )[l-P(0)] 
V D ) Ai-A 2 c[l-P(0)]+z(l- 



c) 



(25) 
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where 



Ai,2 - - 



± 



z(l-c) 



z(l-c) + c[l-V(0)] 



z(l-c) 



4eg[l - V(0)} 



z(l - c) + c[l - V(0)]J z(l - c) + c[l - 7>(0)] 



(26) 



are the eigenvalues of the T matrix in this approximation. From equations ( |25| ) and 
(26) we find that, in the limit p s — > 0, 



Ai. 



and 



1 



1 - c 



I) 



1 -c 



c(l - c) 



py 2 c(Ai - Aj) + gy 2 (c - ^(A^' 1 - A^ 1 ) 
c?(l-c)(Ai-A a ) 



which is the jam-gap distribution in the NaSch model for K^aa = 1 P^j- 

The probabilities like P2(s/c-i, Sfc) which appeared in the intermediate steps of 
the calculations in this subsection are not true 2-cluster probabilities in the usual 
sense of the term The final expression for the distribution Pj g (k) involves only 
quantities which we compute entirely within the framework of the COMF theory. 
Therefore, this whole approach may be called a hybrid approach. We shall compare 
the expression (|2^) for Vj g (k) as well as the corresponding expression derived in the 
next subsection following the true 2-cluster approach with the data obtained from 
computer simulation at the end of this section. 



5.2 The 2-cluster approach 

We now derive the distribution Pj g (k) using the true 2-cluster approximation. Let 
us now allow the state variable s to take one of the four possible values, namely, 
s — — 1 corresponding to an empty site, s — 1 corresponding to a site occupied 
by a vehicle with instantaneous speed V = 1 while s = and 2 correspond to a 
site occupied by, respectively, a non-slow-to-start vehicle with V = and a slow- 
to-start vehicle (also with V = 0). In terms of this new definition of the variable s, 
we can write the distribution of the gaps between the jams in the BJH model with 
'max f 

V jg {k)=M' jg J2 g 2(a | ai)C 2 (gi | s 2 ) ■ ■ ■ C 2 (ffc-i I s k )C 2 (sk | a') (27) 

<r,o-'=0,2 { Si = ±l} 

where Afj g is a normalization coefficient and C2(a|6) are the true conditional 2- 
clustcr probabilities. In terms of the transfer matrix T, whose elements are given 
by the equation 

r[i,i] = c a (i|i), T[i,-i] = c 2 (i|-i), 

T[-l J l] = C a ( = l|l), T[-l J -l]=C 2 ( = l|-l) J 
the equation (|27]) can be written as 

^ g (k)=M; g J2 E C 2 (a\ S )T k ^[ S , S '}C 2 (lW'). (29) 

{cr,cr'=0,2} s,a'=±l 

For convenience we have labeled the elements of T by (1,1), (1,-1)) ("1,1) > ("!,"!) 
instead of (1,1), (1,2), (2,1), (2,2). As before, T fe ~ 1 [z, j] refers to the element of 
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the matrix T fc_1 and are given by the equations (|T7|-pO[) . We compute the 2-cluster 
probabilities as outlined in appendix B and substituting in equation ( p9[ ) get the 
distribution of gaps between successive jams in the BJH model with V max = 1- 
Just as in the case of the jam-size distribution, the agreement of expressions 




Figure 7: The distribution of sizes of gaps between successive jams in the BJH 
model with V max = 1 are shown on a semi-log plot for the same values of c, p and 
p s as in fig. The lines are obtained from the approximate analytical expression 
(29) for Vjg(k) while the discrete data points have been obtained from computer 



simulation. 
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Figure 8: The distributions of the gaps between successive jams in the BJH model 
and the NaSch model, both with V ma x = 1 and p = 0.05, are shown on a semi-log 
plot for three different values of the density of vehicles c. The lines are obtained from 
the exact analytical expressions for Vj g (k) in the NaSch model while the discrete 
data points have been obtained from computer simulation of the BJH model with 
p s = 0.50. 



(25) and ( |29| ) with the corresponding computer simulation data is better when 
the vehicle density is lower (fig. [?]). We have found that, at lower densities, the 
prediction of the hybrid approach (i.e., the equation (p5|)), rather than that of the 
true 2-cluster approximation (i.e., the equation (p9j), is in better agreement with 
the corresponding computer simulation data. However, beyond some intermediate 
density, the results of the true 2-cluster approximation are found to be closer to the 
computer simulation data than the corresponding results obtained in the hybrid 
approach. It is interesting to note that, for the values of the set of parameters 
chosen for fig. |?], there is very little difference between the jam-gap distributions in 
the NaSch and BJH models (see fig. ||). 
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Figure 9: Same as fig. [?], except that V max — 5 and the lines are merely guides to 
the eye. The inset shows the local maximum at £/(c). 

Comparing the curves in fig. [5],|]with the corresponding distributions for V max = 
5, shown in fig. ^, we conclude that the distribution Vj g (k) decays exponentially 
with the gap k between the successive jams for all finite V max , as far as small 
distances are considered. Nevertheless Vj g (k) allows to identify macroscopic free 
flow regimes. In particular, we found for v max — 5, p s — 0.5, and p = 0.05 a local 
maximum of the distance distribution at Lf (c) = L — Lj (c) , where L f (c) denotes 
the average size of the free flow regime at a given density, and Lf(c) the typical 
size of the macroscopic jam. Therefore the jam-gap distribution is well suited to 
identify at least the macroscopic free flow regime of a phase separated high density 
state. 

6 The TH distribution 

We label the position of the detector by j = 0, the site immediately in front of 
it by j = 1, and so on. For the convenience of analytical calculations of the TH 
distribution in the steady state of the BJH model with V max — 1, we assume the 
sequence of steps 2 — 3 — 4 — 5 — 1, just as in the corresponding calculation for the 
NaSch model |l8|,[lj]] the sequence of steps was assumed to be 2 — 3 — 4— 1. 

Suppose, V m (ti) is the probability that the FV takes time ti to reach the de- 
tector moving from its initial position, where it was located when the clock at the 
detector site was set to t = 0, after the LV just left the detector site. Suppose, after 
reaching the detector site, the FV waits there for r — 1\ time steps, either because of 
the presence of another vehicle in front of it or because of its own random braking; 
the probability for this event is denoted by Q(t — t\\t\). Then, the distribution 
Vthir), of the TH r, is given by g§ 

r-l 
ii=l 

or, equivalently, by 

r-l 

Vthij) = V ™{t - t w )Q(t w \r - t w ) (30) 
t w =i 

where t w is the waiting time at the detector site. 

We calculate Vth( T ) f° r Vma X — 1 analytically using the equation (|30"|). In order 
to calculate V m {ti) we consider those spatial configurations at t = from which the 
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FV can reach the detector site within t\ time steps. This implies that one needs to 
consider configurations up to a maximum separation of t\ sites from the detector 
site so that the FV vehicle can reach the detector site even from the farthest point 
with t\ hops. Thus, the configurations of interest are of the form 



(1,-1, 



-11-1) 



where n — 1,2, ...,t\. The underlined — 1 implies that we have to find the condi- 
tional steady-state probability for the given configuration subject to the condition 
that the underlined site (detector site) has just become empty. In the 2-cluster 
approximation, this probability, II(n), is given by 



n(n) = p 2 (i| z i){p 2 (-i|-_i)r- 1 . 

where, the 2-cluster steady-state probabilities are given by 

^2(1|=I) = 3 [1-^0-^(1)] 



(31) 



(32) 



and 



iM-i|=l) = * = 



P9_ 

qg 



(33) 



For all configurations with t\ > n, the vehicle has to stop (t\ — n) times in 
order to reach the detector site after exactly t\ time steps. Equivalently, t\ — 1 time 
steps must elapse in crossing n\ bonds (as the vehicle must certainly hop across the 
bond connecting the detector site with its preceding site at the last time step); the 
number of ways in which t\ — l time steps can be distributed among n — 1 bonds is 

(n-D- ThuS < 



V m {tx) 



£ 

71=1 



tl-1 

n - 1 



n(n)g>* 1_ 



(34) 



as a factor q for every forward step and a factor of p for every halt must be included. 
Substituting j32|) and ( |33"|) into equation (|3l]) and using the corresponding expression 
for n(n) in (p4|) we get 



qc (p 
9 



(35) 



where TVj is the normalization factor. The normalization condition X)t^=i ^m(ii) 
1 leads to 



V m {h 



p (p 

J 



or, equivalently, 



g-p 



l - 



g-p 



It is straightforward to verify that in the limit p s 
hence, from the expression (B6h we find that 



1 _ __ 

d 



0, (g - p)/g 



(36) 

(qy/d) and, 



d 



ti-i 



which is the corresponding expression for the NaSch model pi 
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Next, let us derive an analytical expression for Q(t w \r — t w ). When the FV 
arrives at the detector site exactly r — t w time steps after the departure of the LV, 
the LV can be at any of the sites labeled by 1, 2, ...r — t w + 1. We consider two 
situations separately, namely, (i) the LV is not at site '1' when the FV arrives at the 
detector site, and (ii) the LV is still at site '1' when the FV arrives at the detector 
site. Let us first consider the situation (i). The probability that the LV stays at the 
site T', i.e., it does not move in any of the r — t w time steps = {g) T ~ tm ■ Therefore, 
the probability that the LV is not at site '1' = 1 — {g) T ~ tm ■ If the LV is not at site 
'1' then the probability that the FV halts at the detector site for exactly t w time 
steps is p tm ~ 1 q because it should halt due to random braking for exactly t w — 1 
time steps and move at the last time step. Thus, the contribution to Q(t w \r — t w ) 
coming from all those situations where the LV is not at site T' when the FV reaches 
the detector site is 

[l-igy-^P^-'q (I) 

Let us now calculate the contribution to Q(t w \T—t w ) from those situations which 
are of type (ii), i.e., where the LV is at site T' when the FV reaches the detector 
site. So, the LV will have to move so that the FV is able to leave the detector site 
(and move forward) after t w time steps. Suppose, the LV moves from T' after k 
time steps (1 < k < t w — 1), then the FV will have to stay at the detector site 
for the next t w — 1 — k time steps either because of random braking or because of 
the slow-to-start rule and, then, move forward at the last time step. Note that if 
t w = 1, i.e., if the FV stays at the detector site only for one time step, then the 
site '1' must be empty when the FV reaches the detector site; the corresponding 
contribution to Q(t w \r — t w ) has already been included in (I). The calculation of 
Q(t w \r — t w ) for all r > 2 is quite non-trivial. Suppose, the LV moves from the 
site T' after k time steps and k < t w — 1. Then the FV can stay at the detector 
site during the immediately next time step with probability p s + q s p and for the 
remaining t w — 2 — k time steps with probability p each; in the last step it must 
move forward with probability q. There is, however, yet another possibility. When 
the LV moves from T' after t w — 1 time steps, then the FV has to move forward in 
the last time step; the corresponding probability being qq s . Thus, the contribution 
to all these situations of type (ii) is 



t„-2 



£ (g) k -V-- 2 - k ( Ps + q sP ) + (gY-- 2 q s 

a(gY-- 2 - (3p^- 2 - 



k=l 



g-p 



(ii) 



where a — p s + q s g and (3 — p s + q s p. Thus, combining (I) and (II) we find 

Q(t w \T-t w ) = [l-(gy- 1 ]q for t w = 1 (37) 

and 



Q(t w \r-t w ) = [l-{g) T ~ tw ]p tw ~ x q 



g-p 



+ (g) T ' u "qg 

In the limit p s — » 0, 

Q(r - h\h) -> [1 - (g) tl ] P T -^- l q + (g^gq 



for L„ > 2 



(38) 



(g) 



T-tl-1 



(p) 



T-tl-1 



g-p 



which is, indeed, the corresponding expression for Q(t — ti\t\) in the NaSch model 



15 



Finally, using ( |37| ) and ( p8[ ) in Q30|) and carrying out the summation we get 

\g-p - 1 __ r _ a 



3 (l/ff) - 1 3 

(5 - p)<ip t ~ 2 (t - 1) - gqPp T ^{T - 1) - ff3P T 



g -p 



(- 




\3 





(39) 



where A/jft is the normalization coefficient (which is required as the expression ( |39| 
is not exact). From the expression (|39"|) we also find that in the limit p s — > 0, 



c - y 



i _ < K 



T-l 



qy 



1 _w 

d . 



d — y 

p-"- 1 -q 2 (r -l)p T ~ 2 



_c-y d-y 

which is the known exact expression for Vth(T~) in the NaSch model |H 
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Figure 10: The TH-distribution in the BJH model with V max = 1 are shown on a 
semi- log plot for the same values of c, p and p s as in figs. |l|, ||. The lines are obtained 
from the approximate analytical expression (82) for Vthij) while the discrete data 
points have been obtained from computer simulation. 



The TH-distribution for the BJH model with V max = 1 is well approximated by 
the expression ( |39"| ) at low densities of the vehicles but the approximate expression 
(82) deviates more and more from the computer simulation data with increasing 
density (fig. |l0|). A comparison of the TH distributions in the NaSch and BJH 
models, both for V max = 1 and p — 0.05, in fig. [ll], clearly demonstrates that the 
larger is the density c the wider is the difference between the TH distributions in 
these two models. In fact, at a high density, e.g., c = 0.5, the larger TH are much 
more probable in the BJH model than in the NaSch model. This is basically for 
two reasons. Firstly the slow-to-start rule leads to larger jams compared to the 
NaSch model at the same density (see fig. ||) and secondly the downstream velocity 
of the jams is reduced, which also leads to an increase of the typical waiting time 
for a car in a jam. The computer simulation data for the same densities as in 
fig. [l0|,[ll]; but for Vmax = 5 are shown in fig. [l^. From the definition of TH it is 
clear that, for the NaSch model as well as the BJH model, the TH distributions 
must vanish for r < 2 if V max = 1 but need not do so if V max > 1 pq . Except 
for this difference, at low densities, e.g., c = 0.1 the TH distributions Vthij) for 
Vmax = 1 is qualitatively similar to that for V max = 5; both exhibit a maximum 
and fall exponentially with r beyond the most-probable TH. But, there are some 
qualitative differences in the r-dependence of Vth{ T ) f° r Vmax = 1 and that for 
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Figure 11: The TH distributions in the BJH model and the NaSch model, both with 
V max = 1 and p = 0.05, are shown on a semi-log plot for three different values of the 
density of vehicles c. The lines are obtained from the exact analytical expressions 
for Vth(k) in the NaSch model while the discrete data points have been obtained 
from computer simulation of the BJH model with p s = 0.50. 
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Figure 12: Same as fig. [l0[ except that V max = 5. 

Vmax = 5 in the BJH model. Beyond the most probable TH, the TH distribution 
for Vmax = 1 can be fitted to a exponentially decaying function with a single decay 
rate. On the other hand, for V ma x = 5, there are two regimes of r: in the small 
r regime, which corresponds to the "free-flowing" regions of traffic, 7 > th{ T ) can be 
fitted to an exponentially decaying function (with a fast rate of decay) whereas 
in the large r regime, which corresponds to the jammed region of traffic, V t h( T ) 
can be fitted also to an exponentially-decaying function (with a very slow rate of 
decay). This phenomenon of two different regimes of t, which arises from the phase 
separation in the BJH model, is absent in the NaSch model |Ts|| . 

7 Regime of validity of the analytical expressions 

In the preceding section we have observed that the analytical expressions derived for 
the BJH model with V ma x = 1 are is good agreement with the numerical data from 
our computer simulation if the density of vehicles is not high. However, the higher is 
the vehicle density the stronger are those correlations which have been neglected in 
our calculations within the COMF approximation and the 2-cluster approximation; 
these correlations lead to increasing deviation of the analytical result, from the 
corresponding simulation data, with increasing density of the vehicles. We now 




I c = 0.10 (BJH) o 

c = 0.10 (NaSch) + 
i c = 0.25 (BJH □ . 
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make a quantitative estimation of the regime of validity of our analytical results 
derived in this paper. 

For the NaSch model, the 2-cluster approximation, which gives exact results for 
V max = 1 , is known to be equivalent to the car-oriented mean- field (COMF) theory 
which neglects the correlations of the type 



G„, m =P(n,m) -P{n)V(m) 



(40) 



where V(n,m) is the joint probability of finding a DH m followed immediately 
by a DH n. On the other hand, the COMF does not give exact results for the 
BJH model p8| because of non-vanishing correlations of this type. Therefore, we 
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Figure 13: The correlation function Go,o(c) in the BJH model with V ma x = 1 is 
plotted against c for the same values of the parameters p and p s as in the figs. |l|, 
1110,10, and0. 

have explicitly computed the correlation Go,o as a function of the vehicle density, 
c, for the BJH model with V max = 1 through computer simulation (fig. |l3|). This 
correlation falls exponentially with decreasing density and can be neglected if c is 
not larger than about 0.3. This is consistent with all the results summarized in 
figs. |l|, ||, [|, ||, 0, ||, [l0[ and [ll], where we have observed reasonably good agreement 
between the analytical results and computer simulations for c — 0.10 and 0.25 but 
not for c = 0.50. 



8 Summary and conclusion 

There are, at present, several different CA models, each providing description of 
a range of traffic phenomena for which it has been specifically developed. Here 
we have considered the BJH model |25} | which incorporates in a collective way, 
through the stochastic parameter p s , a specific realistic aspect of individual driving 
habits into the minimal CA model suggested in the pioneering work of Nagel and 
Schreckenberg ]Tl|j . 

In this paper we have investigated in detail the spatio-temporal organization 
of the traffic in the BJH model [^5) by calculating the distributions of DH, TH, 
jam sizes and gaps between jams. For V max — 1, we have computed the DH 
distribution using the COMF approach. We have derived analytical expressions 
for the distributions of jam-sizes and gaps between jams as well as that of the TH 
distribution in the BJH model with V max — 1 by a combination of COMF theory 
and a 2-cluster approximation within the SOMF theory. We have compared these 
approximate results with high quality numerical data obtained from our computer 
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simulations. For V max > 1, we have computed all the above mentioned distributions 
only via extensive computer simulations. 

We have analyzed the data for two different regimes, namely, p s <C p and p s 3> 
p. In the regime p s <C p, there is hardly any significant difference between the 
analytical results and computer simulation data for the BJH model. In the regime 
p s ^> p, we have found, in general, good agreement between analytical results and 
computer simulation data for V max = 1 when the density of vehicles is not too 
high. However, the higher is the density the wider is the difference between our 
analytical results and computer simulation data. This is a consequence of the fact 
that the COMF approximation as well as the 2-cluster SOMF approximation made 
for our analytical calculations neglect some correlation which become non-negligible 
beyond a certain range of vehicle density. In order to directly establish this fact 
we have made a quantitative estimation of the regime of validity of our analytical 
results by computing an appropriate correlation function. 

In order to elucidate the effects of the slow-to-start rule on the spatio-temporal 
organization of the traffic, for each set of values of the parameters V max ,p,c, we 
have compared all our numerical data for the BJH model (p s 0) also with the 
corresponding data for the NaSch model (p s = 0). From this comparison, some 
novel features of the spatio-temporal organization of vehicles in the BJH model 
have emerged. 

In the BJH model for V max = 1 no phase separation takes place and spontaneous 
jamming is possible also in the free-flow regime. In this sense the model is similar 
to the NaSch model. This is also reflected in the quantities we have calculated. For 
Vmax = 5, at sufficiently low densities, i.e., in the free-flow regime, the behavior of 
the BJH model is also very similar to that of the NaSch model in the corresponding 
regime. But, at higher densities, phase-separated states appear in the BJH model 
with V max = 5. Therefore, in this case, the jam-gap distribution is qualitatively 
different from that of the NaSch model, e.g., one finds a local maximum at a macro- 
scopic distance corresponding to the length of the free-flow regime. On the other 
hand, the jam-size distribution is not a good indicator of the appearance of this 
phase separation because the large "jam" , which is formed in this phase-separated 
regime contains 'holes' whereas, according to our definition, only a compact cluster 
of vehicles with V — qualifies as a single jam. 

Measurements on real traffic have shown the existence of metastable states |}6|, 
[57j . These are related to the phase separation |2?]] found here and are an indication 
for the importance of slow-to-start effects in real traffic. Slow-to-start rules lead to 
a reduction of the outflow J out from a jam compared to the maximal possible flow 
Jmax P7fl - Experimentally one finds J m ax/Jout ~ 1-5 |p6| , |3^| . In the BJH model 
for fixed V rnax this ratio is determined by p and p s . p s also determines the velocity 
of jam fronts. Therefore it is possible to tune these parameters in order to get a 
realistic description of jam properties. 
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A Details of the hybrid approach 



We now calculate the 2-cluster-like probabilities p 2 (a, b) involved in the hybrid ap- 
proach to the B JH model with V max = 1 immediately after the randomization stage 
from which one can obtain the conditional probabilities p 2 (a|6) etc. 
•calculation of p 2 (— 1, —1): 

If the configuration of a pair of nearest-neighbor sites is (—1,-1) after the accel- 
eration stage, it does not change during the slow-to-start stage, blockage stage and 
the randomization stage and this observation leads to the relation p 2 (— 1, — 1) = 
P 2 (— 1, —1). Using the 2-cluster-like decomposition, we write 

V(n) = P 2 (+l| l)[P 2 tl| - 1)]- 1 P 2 (^1|1) 

for arbitrary n > 2 and, correspondingly, 

V(n + 1) = P 2 (+l| - l)[P 2 t!| - 1)]"P 2 till) 

which lead to V(n + 1)/V(n) = P 2 ( z l\ - 1) = P 2 (-l, -1)/[P 2 (-1, 1) + P 2 (-l, 0) + 
P 2 (-l,-l)]. Using the result P 2 (-l, 1) + P 2 (-l, 0) + P 2 (-l, -1) = 1 -c (because 
the probability that the left site is unoccupied, irrespective of the state of the right 
site, is 1 — c) we get Vin + 1)/P(n) = P 2 (— 1, — 1)/(1 — c). On the other hand, from 
the equation (||) we get V(n + 1)/V(n) = z. Therefore, 

p 2 (-l,-l) = P 2 (-l,-l) = z(l-c) (A.l) 

•calculation of p 2 (— 1, +1): 

If there there is a vehicle at the leading site (on the right) with speed 1 at the 
end of the randomization stage then it must move forward during the next vehicle 
movement stage. Moreover, the vacancy on the left (following) site of the pair 
can be a part of s string of m successive empty sites where m can be any non- 
zero positive integer. Furthermore, since the DH distribution immediately after the 
randomization stage is identical to that immediately after the acceleration stage (as 
no vehicle movement takes place in between these two stages of the same time step 
of updating) 

oo 

p 2 (-l, 1) = eg J2 [V(m) + P(l)] - cg[l - P(0)] (A.2) 

m— 1 

where we have utilized the normalization condition (jl|). 
•calculation of p 2 (— 1,0): 

The calculation of p 2 (— 1, 0) proceeds exactly in the same manner as in case of 
p 2 (— 1, +1) except for the fact that the leading site (on the right) is now occupied 
by a vehicle with speed which, consequently, will be unable to move during the 
next vehicle movement stage; hence, the expression for p 2 (— 1,0) is obtained from 
that for p 2 (— 1, +1) replacing g in p 2 (— 1, +1) by 1 — g, i.e., 

p 2 (-l,0) = c(l-<7)[l-P(0)]. (A.3) 

•calculation of p 2 (0, —1): 

The configuration 0,-1 can arise immediately at the end of the randomization 
stage in three different ways: 

(i) the vehicle in the following site (on the left) is one which is a slow-to-start vehicle 
which, indeed, slows down in the slow-to-start stage; the corresponding contribution 
to the probability p 2 (0, —1) is cp s V(l). (ii) the vehicle in the following site (on the 
left) is one which is a slow-to-start vehicle which, however does not slow down 
during the slow-to-start stage but slows down at the later stage of randomization; 
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the corresponding contribution to ^2(0, —1) is cq s pV(l). In both the cases (i) and 
(ii) the site immediately in front of the leading site (on the right) must be occupied 
in order that the vehicle at the following site of the pair can qualify as a candidate 
for slow-to-start option, (iii) the vehicle in the following site is not a slow-to-start 
vehicle and it slows down at the randomization stage. Moreover, in this case, the 
empty leading site can be a part of string of m successive empty sites where m 
can be any non-zero positive integer; the corresponding contribution to P2(0, —1) is 
c PEm=i V{m)\ = cp[l — V(0)—V(l)]. Therefore, adding all the three contributions 
we get 

P2 (0, -1) = c\p s V(l) + q s pf(l) +p(l - V(0) - V(l))}. (A.4) 

•calculation of p 2 (0,0): 

The configuration 0, can be obtained immediately after the randomization 
stage if and only if it was +1,+1 immediately after the acceleration stage. The 
facts (i) that the following site is occupied by a vehicle, (ii) that there is no empty 
site in front of the following vehicle, and (iii) that the leading vehicle cannot move 
during the vehicle movement stage (because of its speed being at that stage) lead 
to 

p 2 (0,0) = c(l-g)V(0) (A.5) 

•calculation of P2(0, +1): 

We now apply the same arguments as those in the calculation of ^(0, 0) except 
that now the leading vehicle will move during the next vehicle movement stage; 
therefore, replacing 1 — g in the expression for P2(0, +1) by g we get 

p 2 (0,+l)=cgV(0). (A.6) 

•calculation of p2(T —1): 

The configuration +1,-1 can arise immediately at the end of the randomization 
stage in two different ways: 

(i) the vehicle in the following site (on the left) is one which is a slow-to-start vehicle 
which, however neither slows down during the slow-to-start stage nor slows down 
at the later stage of randomization; the corresponding contribution to p2(+l, —1) 
is cq s qP{l). In this case (i) the site immediately in front of the leading site (on 
the right) must be occupied in order that the vehicle at the following site of the 
pair can qualify as a candidate for slow-to-start option, (ii) the vehicle in the 
following site is neither a slow-to-start vehicle and nor does it slow down at the 
randomization stage. Moreover, in this case, the empty leading site can be a part of 
string of m successive empty sites where m can be any non-zero positive integer; the 
corresponding contribution to P2(+1, —1) is cqJ2m=i T{ m ) = C( z[l — ^(0) — ^(1)]- 
Therefore, adding the two contributions we get 

p 2 (+l, -1) = cq[q s V(l) + (1 - T(0) - P(l))}. (A.7) 

•calculation of p2(+l,0) and p2(+l, +1): Since the speed of a following vehicle 
is certainly reduced to zero at the blocking stage if the leading site is occupied by 
a vehicle, it is impossible to have the configurations +1,0 and +1, +1 immediately 
after the randomization stage; therefore, 

p 2 (+l,0) = (A.8) 

and 

p 2 (+l,+l) = 0. (A.9) 



In the special case p s = the 2-cluster-like probabilities ( A.1)-( [A.9| ) reduce to the 
corresponding expressions for the true 2-cluster probabilities derived earlier Jl7| for 
the NaSch model. 
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B The 2-cluster equations 

Using the symbols fi = p s + pq s and v = qq s we can write the exact equations for 
time evolution of the 2-cluster probabilities C 2 (si,s 2 ) in the BJH model in terms 
of the 2-cluster and higher-cluster probabilities in the steady state as follows. 

C 2 (-l, -1) = C 3 (-l, -1, -1) +pC 3 (0, -1, -1) + pC 3 (l, -1, -1) + mC 3 (2, -1, -1) 

+ qC 4 (-l, -1, 0, -1) + pqC 4 (0, -1, 0, -1) + /igC 4 (2, -1,0, -1) 

+ pgC 4 (l, -1, 0, -1) + i/C 4 (-l, -1, 2, -1) + pi/C 4 (0, -1, 2, -1) 

+ /ii/C 4 (2, -1, 2, -1) + puC 4 (l, -1, 2, -1) + 9 C 4 (-1, -1, 1, -1) 

+ OT C 4 (0,-1,1,-1) + M C 4 (2,-1,1,-1)+ OT C 4 (1, -1,1,-1) (B.l) 

C 2 (0,-1) = K7 2 (0,-l)+,iC 2 (2,-l)+K7 2 (l,-l) (B.2) 

C 2 (2,-l) = gC 2 (0,-l) + i/C 2 (2,-l)- g <7 3 (-l,0,-l)-i/C 3 (-l,2,-l) 

+ qC 3 (0, 0, -1) + i/C 3 (0, 2, -1) + qC 3 (0, 1, -1) + g C 3 (2, 1, -1) 

+ gC 3 (l,l,-l)-«C , 3(0,0,-l)-i/C , 3(0,2,-l) (B.3) 

C 2 (l, -1) = qC 3 (0, -1, -1) + 9 C7 3 (1, -1, -1) + vC 3 (2, -1, -1) + q 2 C 4 (0, -1, 0, -1) 

+ quC 4 {2, -1,0, -1) + g 2 C 4 (l, -1,0, -1) + gi/C 4 (0, -1, 2, -1) 

+ f 2 C 4 (2, -1, 2, -1) + gi/C 4 (l, -1, 2, -1) + 9 2 C 4 (0, -1, 1, -1) 

+ quC i {2,-l,l,-l) + q 2 C i {l, -1,1,-1) (B.4) 

C 2 (-1,0) = pC 4 (-l,-l,0,-l)+p 2 C4(0,-l,0,-l) + W C 4 (2,-l,0,-l) 

+ P 2 G\(l, -1,0, -1) + M C 4 (-1, -1, 2, -1) + W C 4 (0, -1, 2, -1) 

+ /i 2 C 4 (2, -1, 2, -1) + Pf iC A (l, -1, 2, -1) +pC 4 (-l, -1, 1, -1) 

+ p 2 C 4 (0,-l,l,-l)+ W C 4 (2,-l,l,-l)+p 2 C 4 (l, -1,1,-1) (B.5) 

C 2 (0,0) = (B.6) 

C 2 (2,0) - P C 2 (0,-l) + MC 2 (2,-l)+pC 2 (l,-l)-pC 3 (-l,0,-l) 

- ^(-1,2,-1) -pC 3 (-l, 1,-1) (B.7) 

C 2 (1,0) = pqC 4 (0, -1,0, -l)+puC 4 (2, -1,0, -l)+pqC 4 {l, -1,0,-1) 

+ qf iC 4 (0, -1, 2, -1) + ^C 4 (2, -1, 2, -1) + m C 4 {l, -1, 2, -1) 

+ pgC 4 (0,-l,l,-l)+^C 4 (2,-l,l,-l)+ OT C 4 (l, -1,1,-1) (B.8) 

C 2 (-l,2) = C 3 (-l,-l,0)-C 4 (-l,-l,0,-l)+pC 3 (0,-l,0)-pC 4 (0, -1,0,-1) 

+ M C 3 (2, -1, 0) - iiC 4 {2, -1,0, -1) + pC 3 (l, -1, 0) - pC 4 (l, -1,0, -1) 

+ C 3 (-l, -1, 2) - C 4 (-l, -1, 2, -1) + pC 3 (0, -1, 2) - P C 4 (0, -1, 2, -1) 

+ pC 3 (2, -1, 2) - M C 4 (2, -1, 2 - 1) + P C 3 {1, -1, 2) - pC 4 (l, -1, 2, -1) 

+ C 3 (-l, -1, 1) - C 4 (-l, -1, 1, -1) + P C 3 (0, -1, 1) - P C 4 (0, -1, 1, -1) 

+ nC s (2, -1, 1) - M C 4 (2, -1, 1, -1) + pC 3 (l, -1, 1) - pC 4 (l, -1, 1, -1) 

(B.9) 
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C 2 (0,2) = (B.10) 

C 2 (2,2) = 0-2(2, 0) + C a (l,0) + C 3 (2, 2) + C 2 (1, 2) -0-3(2, 0,-1) 

- C 3 (l, 0, -1) - C 3 (2, 2, -1) - C 3 (l, 2, -1) + C 2 (0, 0) 

- C 3 (0, 0, -1) + C 2 (0, 2) - C 3 (0, 2, -1) + C 2 (0, 1) - C 3 (0, 1, -1) 

+ C 2 (2,l)-C 3 (2,l,-1) + C 2 (l,l)-C 3 (l,l,-1) (B.ll) 

C 2 (l, 2) = qC 3 (0, -1, 0) - qC 4 {0, -1, 0, -1) + ^(7 3 (2, -1, 0) - ^C 4 (2, -1, 0, -1) 

+ qC 3 (l, -1, 0) - g C 4 (l, -1, 0, -1) + qC 3 (0, -1, 2) - g C 4 (l, -1, 2, -1) 

+ ^C 3 (2, -1, 2) - ^<7 4 (2, -1, 2, -1) + qC 3 (l, -1, 2) - g(7 4 (l, -1, 2, -1) 

+ <zC 3 (0, -1, 1) - g C 4 (0, -1, 1, -1) + ^C 3 (2, -1, 1) - ^C 4 (2, -1, 1, -1) 

+ gC 3 (l,-l,l)- g C 4 (l, -1,1,-1) (B.12) 

C 2 (-l,l) = qC 2 (0,-l) + uC 2 (2,-l) + qC 2 (l,-l) (B.13) 

C 2 (0,1) = (B.14) 

C 2 (2,l) = (B.15) 

C 2 (l,l) = 0. (B.16) 

In order to close these equations, we now make the approximate decomposition of 
the 3-cluster and 4-cluster probabilities in terms of 2-cluster probabilties using the 
formulae 

n , \ C 2 (x 1 ,x 2 )C 2 (x 2 ,x 3 ) 

C 3 (xi,a; 2 ,ar 3 ) = — — 5 — — (B.17) 

2^ s= -iC 2 (x 2 ,s) 



C 2 (xi , x 2 ) C 2 ( x 2 , x 3 ) C 2 (x 3 , x A ) 



J2l = -iC 2 (s,x 2 ) • EL-iCz^s) 



(B.18) 



We like to remark that the so-called Kolmogorov consistency conditions ]3£| 

2 2 

J2 C 2 (-l,s)= C 2 (s,-l) = l-c (B.19) 



s=-l 



and 



[C 2 (0,s) + C 2 (l,s) + C 2 (2,s)}= [C 2 (s,0) + C 2 ( S ,l) + C 2 (s,2)}=c 

(B.20) 



s=-l 



are useful also for the numerical solution of the equations (B.2-B.16). 



C Relation between COMF and hybrid approaches 

In order to point out the equivalence of the hybrid and COMF approaches, we 
sketch here the derivation of the distribution of the jam sizes in the BJH model 
with V m ax = 1 following the hybrid approach and show that we get the same 
expression as we derived in section 4 following the pure COMF approach. 

Decomposing the distribution Vj S (k) into the 2-cluster-like conditional proba- 
bilities p 2 (&\b) we get 

V ]S (k) = ^sP 2 (-A\0)[ P2 (0\0)} k ^P2(0\ - 1) +P2(zl|0)b2(0|0)] fc -V 2 (0|l) (C.l) 

where Afj S is a normalization coefficient. Using the expressions for the 2-clustcr-likc 
probabilities derived in the appendix A and normalizing the distribution to unity 
we recover ]35| the expression ( |T2| ) of Vj S (k) in terms of 'P(O). 
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